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Abstract 

The role of .S-mixing in the quantum tunneling of the magnetization in nanomagnets has been 
investigated. We show that the effect on the tunneling frequency is huge and that the discrepancy 
(more than 3 orders of magnitude in the tunneling frequency) between spectroscopic and relaxation 
measurements in Fes can be resolved if 5-mixing is taken into account. 
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Molecular nanomagnets (MNMs) are molecules containing transition-metal ions 

whose spins are so strongly exchange-coupled that at low temperature each molecule behaves 
like a single-domain particle with fixed total spin. One of the most interesting phenomena 
displayed by these systems is quantum tunneling of the direction of the total spin through 
energy barriers^, U, 0]. The measured step-like magnetization curves of Mn 12 and Fe§ 
provided macroscopic evidence of relaxation through quantum tunneling. The latter is 
revealed by resonances observed in the relaxation rate at specific values of the external 
magnetic field Bac, at which energy levels on opposite sides of the anisotropy barrier are 
nearly degenerate and anticrossings (ACs) in the field dependence of the energies occur. The 
relaxation rate depends crucially on the value of the so-called tunnel splitting A, i.e. the 
gap at Bac between the quasi-degenerate states. In particular, at very low temperature T 
and at short times t the magnetization relaxes as 1 — VTi where T oc A 2 [6]. For Feg, A 
was extracted in Ref. by measuring with a microsquid apparatus the magnetization steps 
induced by sweeping a longitudinal (i.e. parallel to the easy axis) applied field B z across 
Bac- The size of these steps was linked to the tunnel splitting through the Landau-Zener 
formula 

n = 1 - e - A2 / A , (1) 

where II is the tunneling probability at a level anticrossing, and A is proportional to the 
field sweeping rate. When the experiment is performed in a static transverse field B y , A is 
found to display oscillations as a function of B y , which in a semiclassical approach reflect 
the destructive interference of tunneling pathways. 

Eq. (1) had been deduced by neglecting decoherence sources such as hyperfine and dipolar 
fields. Nevertheless, it remains valid if the sweeping rate is as fast as that actually used in 
the experiments (3]. In addition, the model proposed in Ref. [a] shows that the incoherent 
Zener tunneling can be described by Eq. (1) with A renormalized by a factor 

A striking circumstance is that the measured value of A seems completely incompatible 
with the value calculated by using the Hamiltonian determined by inelastic neutron scatter- 
ing: (INS) 0,0, optical spectroscopy llj and electron paramagnetic resonance |l2j|. Indeed, 



the measured zero- field gap A(B y = 0) between the two lowest levels is near 10 7 K0], while 
the value calculated (4.44 ■ 10" 11 K) is more than 3 orders of magnitude smaller. This huge 
discrepancy seriously hinders any attempt to reach a satisfactory theoretical modelling of 
the quantum tunneling of the magnetization, where the square of A plays a fundamental 



ected quantum fluctuations of 



14J) affect the tunnel splitting 



role. The purpose of this work is to show that commonly nee 
the magnitude of the total spin of the molecule (S-mixing |l3j 
of Fe 8 hugely, and allow the above-mentioned discrepancy to be solved. Since Fe 8 displays 
a relatively small degree of .S-mixing, we expect the tunnel splitting of many nanomagnets 
to be influenced even more heavily than in Fe§ by such fluctuations. 

Each Fes molecule can be described by the following spin Hamiltonian: 

H = J2 Jii* ■ Si + E E b l(i) l( s i)+ 

i>3 i k,g 

E s i ■ D *i • s j - 9Vb E B ' Si ' ( 2 ) 

i>j i 

where Sj are spin operators of the i th Fe 3+ ion in the molecule (s, = 5/2). The first term is 

the isotropic Heisenberg exchange interaction. The second term describes the local crystal- 

rn 

fields (CFs), with 0^(s ; ) Stevens operator equivalents for the z-th ion|l6j and b q k {i) CF 
parameters. Here k = 2 or 4 (larger values of the rank k are forbidden for (i-electrons 16j). 
and q = —k, k. The third term represents the dipolar anisotropic intra-cluster spin-spin 
interactions. The last term is the Zeeman coupling with an external field B. The exchange 



constants Jij used in this work are those determined from susceptibility 

While the Heisenberg term is rotationally invariant and therefore conserves the length 
|S| of the total spin S = ^Sj, the anisotropic terms do not conserve this observable. 
Nevertheless, since the Heisenberg contribution is usually largely dominant, |S| is nearly 
conserved, and the energy spectrum of H consists of a series of level multiplets with an 
almost definite value of |S|. By neglecting the mixing between different S'-multiplets (i.e. S- 
mixing), the Hamiltonian Eq. (2) can be projected onto each S'-multiplet (strong-exchange 
limit): 

H sub = B° 2 0° 2 (S) + Bl0 2 2 {S) + £? 4 O 4 °(S) 

+BjOj(S) + BlO\{S) - gf i B B ■ S, (3) 

where S is a vector spin operator with S equal to the total spin of the S'-multiplet Q- Th e 
parameters B® are calculated from b 9 k (i) and D^- by CF and dipolar projection coefficients. 
This approach, applied to the S — 10 ground manifold of Fes, has allowed to interpret INS 
data very satisfactorily by assuming B% = —9.75 x 10~ 2 K, B\ = —4.66 x 10~ 2 K,B 



1.0 x lO-^Ef = 1.2 x 10- 7 K,Bj = 8.6 x lO -6 ^ Q, Q- Ver y similar parameter values 
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are obtained from optical spectroscopv|ll|| and electron paramagnetic resonance |l2j. In 
particular Bf, which has the greatest effect on A, is the same. In order to reproduce the 
measured magnitude and oscillations of A with Eq. (3), values of B^ incompatible with 
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neutron results had to be assumed |5l. I8L Il9j|. In particular, B\ was one order of magnitude 
larger and its sign was reversed. In Fig. 1 we show calculations at T = 9.6/T of the INS 
spectrum of a Fes powder with the experimental resolution of Ref. jjj. The parameters used in 
Refs.flflQ do not reproduce the INS spectra satisfactorily, neither in the higher-energy 
part measured in Ref.|9| nor in the lower-energy part measured with the high resolution 
experiment of Ref . [lf| . If the strong-exchange-limit Hamiltonian (3) is used, there is no 
way to reproduce with a unique set of parameters the magnitude and behavior of A and 
the spectroscopic results. In the following we show how this discrepancy can be removed if 
^-mixing is taken into account. Indeed, although in Fes S'-mixing is a small perturbation 
(e.g., it produces negligible changes in calculated spectroscopic quantities), its effect on A 
is very large because it provides much more efficient tunneling channels. 

In order to evaluate S'-mixing effects, we followed the method developed in jl3jj], in which 
S-mixing is included up to the second order in the anisotropy by a unitary transformation 
applied to the complete Hamiltonian (2). The system can be still described as an effective 
spin S = 10, provided the spin-Hamiltonian (3) is properly modified : the parameters of the 
Stevens operators are renormalized, and new higher rank (K > 4) terms are added. These 
latter are forbidden for <i-electrons in the strong-exchange limit. The advantage of using 
this method with respect to large-scale numerical diagonalization (e.g. using the Lanczos 
algorithm) is twofold: first of all, it allows calculation times to be reduced drastically. In 
fact, the time-consuming part of calculation (i.e. computing H and T in Eq. (9)) does 
not depend on the specific set of local CF parameters, and therefore has to be performed 
only once. The second advantage is that the simple and physically transparent single-spin 
formalism of the strong-exchange limit is recovered. 

Using as basis vectors the eigenvectors \aSM) of the isotropic exchange Ho, the full 
Hamiltonian matrix H Eq. (2) can be written as the sum of three terms 

H = H + H 1 + H 2 (4) 

where Hi + H 2 represents the anisotropic interactions. H\ has nonzero elements only within 
the S multiplets, while H 2 joins states with different aS and is the term responsible for the 
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mixing. H 2 is neglected in the strong-exchange limit Eq. (3). 
The perturbational procedure 



13l . consists in performing a unitary transformation 
on H such that the off-diagonal (in aS) blocks of the transformed Hamiltonian H' are zero 
up to second order in the anisotropy. Hence, in the new basis, states belonging to different 
multiplets are uncoupled and the system can be described as an effective spin multiplet, like 
in the strong-exchange limit. The matrix elements of H' inside the ground multiplet S = 10 
are given by 

(SM\H'\SM f ) = E q 5 m ,m> + (aSM\H 1 \aSM') 

x ^ (aSM\H 2 \a"S"M")(a"S"M"\H 2 \aSM') 

~ 1^ p F ^ 

a"S"M" 

where Eq is the lowest eigenvalue of Hq and \aSM) are the corresponding eigenvectors. 
\ot"S"M") are excited eigenvectors of Hq with energy £? a"S"- 

The second term in Eq. (5) coincides with the strong-exchange Hamiltonian (3), while 
the last term represents mixing corrections. By exploiting the Wigner-Eckart theorem 



the latter can be written in general as 



Y.Bk°k (6) 

K,Q 

with K < 8 and even, and —K < Q < K. Hence 

where C% = B% + B% (B% = for K > 4). Therefore, on the one hand S'-mixing introduces 
in the effective Hamiltonian new terms forbidden in the strong-exchange limit, and on the 
other hand it renormalizes the coefficients of the other terms (with K < 4). This implies 
that the CF parameters determined by INS are to be regarded as C^s rather than as the 
S^sofEq. (3). 

The parameters are given by linear combinations of products of reduced matrix 
elements. For example, 

Bt = E & 2«^(J)2(^,J) + E & 4«TW (8) 

S+2 

s(<,i)= E 



with 



5+2 

r S~S" 
c 66 



y (aS\\T^ma"S")(aS\\T^(j)\\a"S") (Qq) 

„ E 0a "S" — Eq 



and m = -VT£j lk £ 



C 66 



y (^S\\T^U\jkW'S^ ^ {%) 

„ E 0a "S" ~ E 

a" 

i, j and /c label magnetic ions, T^ K \i) and T^(ll\jk) are the tensor operators describing the 
local CF and dipole-dipole interactions [20]. J™ k are defined in terms of the elements of Djj 
(Ea. (2)) in Ref. j^. The Cg 6 ~ 5 coefficients are defined according to the theory developed 
in jl3l] . Expressions similar to (8), (9a) and (9b) are obtained for the other parameters B®, 
which are all expressed as polynomials of 2 nd -order in the b q k (i). Hence S- mixing gives rise 
to highly efficient and otherwise forbidden tunneling channels, by generating new high-rank 
anisotropy terms. 

In order to assess the impact of these terms on the tunnel-splitting we applied our theory 
quantitatively. While the dipole-dipole interaction (appearing, e.g., in (9b) through J ?l ) can 
be computed by the point-dipole approximation [2JJ, the local CF parameters b q k {i) cannot 
be determined ab initio reliably. Therefore, by numerically inverting the 2 nd -order functions 
C®({b q k (i)}) , we determine the possible sets {b q k (i)}f consistent with INS, i.e. such that the 
values C^({b1(i)} f) f° r K <A coincide with those determined by INS (within experimental 
error bars). Even by neglecting all b q k (i) with q ^ 0,2,4 (i.e. those not contributing to Hi) 
and by enforcing on the {b k (i)} the approximate D2 molecular symmetry of Fes [22(, there 
are still more unknown parameters than constraints, and we find therefore that there are 
infinitely many sets compatible with INS. For these sets, the distribution of the calculated 
values of the tunnel splitting A is shown in Fig. 2a by a histogram of the log-increments 
log(A/A ). Here the local CF parameters b q k {i) vary on a grid bounded by \b\{i)\ < 8 K, 
|&4(i) I < 0.4 K. This choice is based on two considerations: i) the experimental values of G\ 
and Cf set only a lower bound I&2WI ^ ^-K. Local second-order parameters of the order of 
few K are reasonable in case of Fe 3+ in a low symmetry environment j^. ii) Typical ratios 
of fourth- to second-order CF parameters range (in modulus) between 0.01 and 0.1. Fig. 
2a shows that S"-mixing plays a crucial role since typically it enhances A by several orders 
of magnitude. Fig. 2b reports the result of the same calculation when the grid bounds are 
restricted to \b\{i)\ < 4 K, \b\{i)\ < 0.04 K. Even with this more restrictive choice the effect 
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of S- mixing remains huge. The measured value of A is indicated by an arrow and falls well 
inside the distribution. 

We have shown that S"-mixing can remove the discrepancy between the value of the zero- 
field gap A(B Z = 0, By = 0) measured by relaxation experiments and that calculated from 
the spectroscopic results. Now it remains to prove that also the measured oscillations of 
A(B Z = 0, By), and the behavior of the AC gaps A ex {B z = Bac, B y ) between excited states 
are reproducible. The aim of this work is not to perform a best-fit of the observed oscillations 
of the tunnel-splittings, but to prove that S"-mixing eliminates the inconsistency between 
spectroscopic and Landau-Zener measurements. Therefore, we limit Eq.(6) to values of 
K < 6 in order to find among the infinite possible parameter sets consistent with INS, 
one involving as few high-rank terms as possible, and reproducing the AC gaps behavior 
satisfactorily . The new terms in Eq.(7) (forbidden in the strong-exchange limit) with 
significative influence on A are then C| and C|. 

As a first step, we fixed for K < 4 to the values determined by neutron spectroscopy 
(reported below Eq. (3)) and we chose values of C| and C| which reproduce the behavior of 
A satisfactorily. With C\ ~ -1.8 x 1(T 7 K and C\ ~ -1.15 x 10~ 7 K, A calculated at zero 
applied field is ~ 1.1 x 1CT 7 K, to be compared with the value ~ 0.4 x 1CT 10 K obtained 
when Cq = and Cg = 0. Moreover, the measured oscillations of A as a function of the 
transverse field B y are well reproduced, as well as the AC gaps involving the excited states 
| — 10) and 1 10 — n) (n = 1,2) (see Fig. 3). The nonzero value of experimental oscillations 
at minima may arise from experimental artifacts, (e.g., crystal mosaicity). The behavior of 
A as a function of the transverse field modulus for nonzero azimuthal angles <fi between the 
applied field and the y-axis, is also in good agreement with measurements^. 

As a second step, we checked that the addition of these 6th-order terms does not affect 
the INS cross-section significantly. In fact, the recalculated cross-section is indistinguishable 
from that calculated in |2j and reported in Fig.l, apart from an irrelevant shift (by ~ 20 
fieV) of the shoulder at 0.16 meV. 

Values of C| and Cg of the order and sign of those given above are realistic in Fes- 
For example, the insets in Figs. 2a and 2b show the distribution P(Cq) of values of Cf 
calculated on the same grids as described above. We stress that our choice of high-rank 
parameters is merely the simplest possibility. There are many different sets involving also 
the other high-rank terms which would be consistent with experimental data. Although a 
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unique determination of the Cj^s is not possible, the important point is that the addition 
of high-rank terms, which is allowed only if S'-mixing is considered, is essential to describe 
consistently relaxation and spectroscopic data. 

In conclusion, we have shown that the discrepancy (more than 3 orders of magnitude 
in the tunnel splitting) between spectroscopy and relaxation measurements in Fe§ can be 
resolved if S-mixing is taken into account. Even a small degree of S-mixing has huge influence 
in the tunneling dynamics since it opens highly efficient tunnel channels through otherwise 
forbidden high-rank anisotropy terms. The degree of S'-mixing is strongly influenced by the 
topology of the molecule. Therefore in addition to the height of the anisotropy barrier, also 
the cluster topology must be taken into account in designing new nanomagnets. 
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FIGURE CAPTIONS: 

Fig.l: Calculated INS intensity for a Fes powder with the spin Hamiltonian Eq. (3) and 
variousparameter sets. The energy resolution has been fixed to the experimental value of 
19/j.evj9(. The parameters used in Kef.! 10 are close to those used in Ref. Q]and yield almost 
the same intensity curve. The inset shows a schematic view of Fe§. 

Fig.2: Calculated distribution of values of the tunnel splitting A(B Z = 0, B y = 0) nor- 
malized to the value Ao = 4.44 ■ 10~ n K obtained without S-mixing (with the parameters 
obtained from INS). The local CF parameters b q k (i) vary on grids (different in (a) and (b)) 
defined in the text. Arrows indicate the measured ratios. Insets show the distribution of 
values of C| on the same grids. 

Fig.3: Top: measured tunnel splitting as function of an applied transverse magnetic field 
By with B z — (n — 0), and AC gaps involving the excited states | — 10) and |10 — n) 
(n = 1,2) with B z = B^c ~ ^-0.22 T (n = 1, 2) p. Bottom: the same quantities calculated 
with the Hamiltonian (7) and the parameters given in the text. 
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